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Abstract 

We study a 2-body problem given by the sum of the Newtonian po- 
tential and an anisotropic perturbation that is a homogeneous function 
of degree —(3, (3 > 2. For (3 > 2, the sets of initial conditions leading 
to collisions/ejections and the one leading to escapes/captures have 
positive measure. For (3 > 2 and /? 7^ 3, the flow on the zero-energy 
manifold is chaotic. For (3 = 2, a case we prove integrable, the infin- 
ity manifold of the zero-energy level is a disconnected set, which has 
heteroclinic connections with the collision manifold. 



PACS(2003): 05.45.-a, 45.50. Jf, 45.50.Pk 
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I. Introduction 



In the past three centuries, celestial mechanics has stimulated the develop- 
ment of many branches of mathematics (see e.g. Refs. 11,13). This trend 
continues, and even its most basic question, known as the 2-body problem or 
the Kepler problem, still attracts the vivid interest of mathematicians and 
physicists, both in its classical form (see Ref. 1) and in newer versions. 

Among the latter are the problems raised by quasihomogeneous potentials, 
given by the sum of homogeneous functions, and problems set in anisotropic 
spaces, for which the interaction law is different in each direction of the space. 
For quasihomogeneous problems the terminology and the first qualitative 
results were introduced in the mid 1990s;^'^^'^^ this potential unifies several 
dynamical laws, including those of Newton, Coulomb, Manev, Schwarzschild, 
Lennard- Jones, Birkhoff and others. The anisotropic case is more related to 
physics, and was initiated by Martin Gutzwiller in the 1970s^'^'^^ for the 
quantization of classical ergodic systems. Among Gutzwiller's goals was also 
that of finding connections between classical and quantum mechanics. A 
combination of the quasihomogeneous and anisotropic aspects shows up in 
the anisotropic Manev problem, whose dynamics contains classical, quantum 
and relativistic features (see Refs. 6,10,14,15). 

In the present paper we consider a new version of the Kepler problem, 
which combines two of the above characteristics: isotropy and anisotropy. 
The potential (see formula ([3])) is the sum of the classical Keplerian potential 
and an anisotropic perturbation, the latter being an homogeneous function of 
degree — /3, (3 >2. This is the first analysis of a quasihomogeneous potential 
that mixes isotropic and anisotropic components. For previously studied 
problems, all terms have been either isotropic or anisotropic. As we will see, 
this case has some surprising dynamical properties, often very different from 
the ones that characterize potentials whose terms are not mixed. 

Such mixed potentials can be used to understand the dynamics of satel- 
lites around oblate planets or the motion of stars around black holes. Here, 
however, we are not interested in applications. Our endeavours are restricted 
to mathematical results. 

In Sec. II we introduce the basic notations, the equations of motion, 
and put into the evidence the symmetries of the problem. In Sec. Ill we 
begin the study of the case /9 > 2, define the collision manifold, which is 
an essential qualitative tool, and perform a geometric study of the flow in 
the neighbourhood of collisions. We classify all collision-ejection orbits and 
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prove that the set of initial conditions leading to them has positive measure. 
We achieve this while studying the flow on and near the collision manifold 
in terms of McGehee-type coordinates. 

In Sec. IV we investigate the existence of heteroclinic orbits on the col- 
lision manifold for potentials with fi > 1 and (3 = 2 + or (3 = 2 + 
The main result of this section is that for an open and dense set of //-values, 
saddle-saddle connections do not exist on the collision manifold. Section V 
deals with capture and escape orbits in the zero-energy case. We show that 
the infinity manifold has two circles of normally hyperbolic equilibria, one 
attractive and one repelling. This proves the existence of infinitely many 
capture and escape orbits. 

In Sec. VI we consider the case /3 = 2, which we show to be integrable. 
Apparently this is quite surprising since the anisotropic Manev problem, 
which resembles this case except for the anisotropy of the Newtonian term, 
is nonintegrable (see Refs. 3,15). But as we will show, the surprise element 
vanishes once we look at the problem in the larger context of the Hamilton- 
Jacobi theory. In Sec. VII we study the flow on and near the collision mani- 
fold and see that its qualitative behaviour is similar to that of the anisotropic 
Manev problem. 

In Sec. VIII we prove the existence of heteroclinic orbits connecting the 
collision and infinity manifolds in the zero-energy case for (3 = 2. The in- 
finity manifold is formed by two disjoint circles of equilibria, which seems 
very unusual. In fact this is the first case we have ever encountered in which 
the infinity (or collision) manifold is a disconnected set. A complete under- 
standing of the topological structure of the phase space when all parameters 
are varied will probably shed more light in this matter. But a study in this 
direction, though interesting and worthy of further investigations, is beyond 
our present scope. 

In Sec. IX we consider a perturbative approach of the problem. The 
perturbation function of the Hamiltonian is a homogeneous function of degree 
—(3 with (3 > 3/2. We end the paper with Section X, in which we apply the 
Melnikov method to show that for every (3 ^ 2,3, the flow on the zero-energy 
manifold is chaotic. This also proves the nonintegrability of the problem for 
(3^2,3. 
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II. Equations of Motion and Symmetries 

Consider the Hamiltonian 

Hpiq,p) = ^p2 + f/^(q), 
where q = {x,y) and p = (px.Py)- The equations of motion are 



(1) 



q= p 

p = -Vt/^(q), 
where Up is a. potential of the form 



(2) 



1 



(3) 



with the constants beta > 2, fi > 1 and 6 > 0. The symmetries of ([2]) 
are given by the following analytic diffeomorphisms in the extended phase 
space: 



Id 

So 
Si 

S2 

S3 
S, 

S5 
S, 



{x,y,Px,Py,t) 
{x,y,Px,Py,t) 
{x,y,px,Py,t) 
{x,y,Px,Py,t) 
{x,y,Px,Py,t) 
{x,y,Px,Py,t) 
{x,y,Px,Py,t) 



{x,y,px,Py,t), 
{x,y, -px, -Py, -t), 
(x, -y, -px,Py, -t), 
{-x,y,px, -Py, -t), 
{-x, -y, -px, -Py,t), 

-Px,Py,t), 

{x, y, Px^ Py^ t^^ 
-y,Px,Py, -t), 



(4) 



where Id is the identity. These diffeomorphisms form a group that is isomor- 
phic to Z2 X Z2 X Z2 (see Refs. 14,22 for more details). Invariance under 
these symmetries implies that if 7(t) is a solution of ([2]), then also 6*4 (7 (t)) 
is solution for i G {0, 1, 2, 3, 4, 5, 6}. 



III. The Collision Manifold for > 2 

We will further express the equations of motion in McGehee-type coordinates, ^'^^ 
which are suitable for understanding the motion near collision. The trans- 
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formations are given step by step as follows. Take 



r 


2 , 2 

= X +y , 




e 


y 

= arctan— , 






X 




V 


= rr = {xpx 4- 


-yPy): 


u 


= r^9 = {xpy ■ 


- ypx. 




13-2 _ 


<3-2 _ 



(5) 



rescale v and u by 



and then rescale the time variable using the transformation 

dt 
dr 

The equations of motion take the form 
r' = rv, 



^/3/2+l_ 



2 A/3/2 ' 

B' = u, (6) 
p-2 6/5(u- l)sin2^ 

U = ——UV+ -3^5 , 

where A = /i cos^ 6' + sin^ 6' and the prime denotes differentiation with respect 
to the new independent time variable r. For simplicity, we keep the same 
notation for the new dependent variables. 

In these new coordinates, the energy integral Hjs = h (see equation ([1])), 
takes the form 

«2 + ^2_2r/3-i- JL = 2/ir^. (7) 
We define the collision manifold (see Figure 1) as 

2h 

C = {{r,v,9,u)\ r = 0, u' + v' = ^}. (8) 

Notice that C is homeomorphic to a torus. The flow on this manifold is given 
by the system 

e' = u, (9) 

, p-2 bp{n-l)sm2e 
" = 2A^ • 




Al/2 



Figure 1: The collision manifold C 

We can now prove the following results. 

Proposition 0.1. All the equilibrium points of system ^ lie on the collision 
manifold C and are given by 



2b 

r = 0, w = ±^/— , ^ = 0,7r/2,7r,37r/2, u = 0. 

Proof. It is obvious that the above points are equilibria for the flow defined 
by (jni). To check that there are no equilibria outside the collision manifold, 
from the first equation in ([H]) with r 7^ we see that v = and from the 
third equation in (Q, that u = 0. Substituting these values in ([7]), we obtain 

Solving the equation = in ([6]), we are led to 

r^-^ + 2hr^ = ^-^^. (11) 

But ffTOj) and ffTTl) have no common solutions, therefore there are no equilibria 
outside the collision manifold. □ 
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Proposition 0.2. The flow on the collision manifold C is gradient-like (i.e. 
increasing) with respect to the {—v)- coordinate. 

Proof. Since /? > 2, we see from the first equation in (Q that t;' < except 
at equilibria, therefore the flow on C increases with respect to —v, so is 
gradient-like relative to it. □ 

To match the sign of v and the value of 6, we denote the equilibria on C 
by , respectively. Observe that A(0) = A(7r) = /i and 

A(-7r/2) = A(37r/2) = 1. With this notation, we can describe the following 
properties of the flow. 

Theorem 0.1. On the collision manifold C , the equilibria and are 
saddles, A^,„ and At ,„ are sources and A~ ,„ and A~ ,„ are sinks. Out- 
side C , the equilibria Aq , A^ , A'^j^ and A^^^j^ have a local one- dimensional 

unstable analytic manifold, whereas Aq , A~ , A^j^ and A'^^i^ have a local 
one- dimensional stable analytic manifold. All these equilibrium points are 
hyperbolic. 

Proof. Consider the function 

26 

F(r, V, e, u)=u^ + v'^ - r^~^ - — — - Ihr^ = 0. 

According to equation ([7]), the surface of constant energy Mh defined by the 
equation 

F{r,e,v,u) = 

is a 3-dimensional manifold. At every point B of M^, the tangent space is 
given by 

TbF = {(r, V, e, u)\VF{B) ■ (r, v, 9, u) = 0}. 

At any equilibrium point A, the tangent space is defined by 

TAF = {{r,e,v,u)\v = 0}. 

A straightforward computation shows that at the equilibria A^ and A^ 
the linearized system corresponding to ([6]) has the matrix 




(/5-2)±(y^) 
1 

±^(V^) 
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therefore the hnear part of the vector field restricted to T^J ^ is given by 

/ 

L = 



± 



26 „ 
772^ 



\ 



u 



\ 



As a basis for the tangent space T^^l~, we take the vectors 



/ 1 \ 









1 






V 1 y 



:i2) 



The representation of the linear part L relative to this basis is given by the 
matrix 



J* 



2b 



v 











fJ+2 





1 



■ (/3-2), 



2b ' 

,13/2, 



\ 
I 



The characteristic polynomial shows that all eigenvalues are real and that 
the equilibrium is a saddle in each case. 

At the equilibria ^^/2 ^^"^ '^tn/2^ same linearized system has the 
matrix 



±V26 

{(3-2)±{V2b) 

1 

-b/3{fi - 1) ±i^{V2b) J 



Using the vectors given in (fT2l) as a basis for the tangent space 



linear part is given by the matrix 



J* 



7r/2,37r/2 ' 



the 





1 

■ (/3-2)/ 



2b) 



8 



The eigenvalues at A'^^^ and at A^^^^ are 



2b, ^^^(v^) + J^P-2)2-8/?(/.-l)], 



and 



^^^-/26)-.r[(/5-2)2-8/5(/i-l)]. 



2 ' ' V 2' 

This means that all of them are positive or have positive real part. At 
and ^3^/2 eigenvalues are all negative or have negative real part. □ 

Corollary 0.1. The set of initial conditions leading to collisions or ejections 
has positive measure. 

Proof. The equilibrium points A~^^ (^^72) ^37r/2 (^3^/2) sinks (sources) 
for the global flow, therefore their basin of attraction (repulsion) is a 3- 
dimensional set of collision (ejection) orbits. □ 

Remark 0.1. If < < , all the eigenvalues are real and positive. 

If /i > ^ , there are two complex eigenvalues with positive real part, so 
some orbits have the spiraling property, i.e. engage into an infinite spin. For 
example, to have spiraling orbits in the case /3 = 3, it is necessary that 
/i > ||. Therefore when /i is sufficiently close to 1, no spiraling orbits exist. 

Remark 0.2. The inequalities in Remark 3.1 do not depend on the param- 
eter b. 

Remark 0.3. For /i close enough to 1, though the set of collision orbits has 
positive measure, there are no spiralling orbits. 



IV. Saddle-Saddle Connections on C. 

Using ideas similar to those found in Ref. 7, we will further study the exis- 
tence of saddle-saddle connections on C for /i > 1 and (3 of the form 

[3 = 2 + -^— and /? = 2 



2k "^^ l + k' 

k integer, k ^ —1. The reason for choosing these values of (3 will be clarified 
below. In the proof we restrict ourselves to the cases /3 = 3, 4. This is because 
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the method requires the computation of an integral for each value of k, and 
each integral has to be computed separately. But the principle is the same 
for every such k. 

To show that there are no heteroclinic connections is sufficient to prove 
that the stable and unstable manifolds of corresponding fixed points miss 
each other. We will show that this holds true for most values of /i > 1. 

It is now convenient to introduce different coordinates on C. Since C 
is homeomorphic to a torus, we can describe the flow using angle variables. 
With the transformations 



.2b 

«=X^sm^, 

r- 13) 

^=^cosV', 

we can rewrite the flow on the equations of motion ([9]) on C as 



ijj =^—^^smiP + -^^sm2e cos ifj. 

The equilibrium points in the new variables {0,tp) are ^^7^/2 ~ (i^/2,0), 
= (±V2,7r), Aq- = (0,0), = (0,7r), A- = (vr, 0) and At = (71,71). 
Notice that if /i = 1 (the isotropic case), the collision manifold is a torus 
for which the upper and lower circles and C~ consist of fixed points. The 
equations take the form 



e' =V2bsmip 

13-2 ^ (15) 
if)' = — ^ — y2hsm.%lj. 

It ie easy to see that in this case there are heteroclinic orbits that connect 
the critical points Aq to A~ and A^. to Aq when 

2 



l + 2k'' 



k integer, 7^ — 1, and connect the critical points A^t^ to A^ and A^ to A_ 
when 
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(a) 



(b) 



Figure 2: The collision manifold and the heteroclinic connections for (a) 
/i = 1, = 3 and for (b) /i = 1, /? = 4. 



k integer, k ^ —1. 

Figure 2(a) depicts the collision manifold for /i = 1 and (3 = 3 and the 
heteroclinic connections on it, while Figure 2(b) does the same for /i = 1 and 



In the following we will show that if/i — l = e>0 and small, such saddle- 
saddle connections are broken and the same result holds for an open dense 
set of /i — 1 = e > 0. 

Theorem 0.2. For j3 = 3 and for an open and dense set of real numbers 
yU > 1, the unstable manifolds at AZ^^ and A~ miss the stable manifolds at 
and A~^^. For (5 = 4: and for an open and dense set of real numbers fi > 1, 
the unstable manifolds at Aq and A~^^ miss the stable manifolds at A~ and 
A-,. 

Proof. Eliminating the time from equations (1140 and using the fact that 



(3 = A. 




(16) 
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we have 



-de=d-r[Am) + [-Y- ) = ^' ^''^ 

where e = fj, — 1 and A = 1 + e cos^ 9. 

First consider /? = 3, and the unstable manifolds W^{—7r, 0) = W^(vr, 0). 
When e = 0, W|(vr, vr) matches iy3*'(— tt, 0). Consider the branch of W^{—7r, 0) 
which contains (0,7r/2). This curve lies along the line 

-2tp + e = -7r. (18) 

When e varies, this branch of the unstable manifold W^{—7t, 0) varies smoothly 
on C. Let C^{9, e) denote the t/^— coordinate of this curve, with (^{—n, e) = 0. 

Now let /3 = 4 and W^{-n,0) = W^{n,0). When e = 0, WI{-tt,0) 
matches PV^ (0, tt). Consider the branch ofW^{—7i, 0) that contains {—n/2, 7r/2). 
This curve lies along the line 

-2^ + 29 = -2n. (19) 

As e varies, this branch of W^{—7i,0) varies smoothly on C. Let C^{9,e) 
denote the coordinate of this curve, with ('^{—ir, e) = 0. Now we need the 
following result, which we will prove at the end of this demonstration. 

Lemma 0.1. Wtth the above notations, £C^(0, 7i/2) = 3/4n > and £C''(-7r/2, 
Tv/2 > 0. 

From this lemma follows that (^{0, e) > and C^(0, e) > for e > small. 
Thus, it is easy to show that 1^1+^(0, C'(0, e)) > 0, where / = 3,4. Equations 
f ll7p are reversed by the transformation 

{9,tP)^{-9,7r-^). (20) 

If /3 = 3, the unstable manifold through {—it, 0) is mapped onto the stable 
manifold through (vr, vr). Hence the stable manifold intersects the line 9 = 
at some point {0,ipo) such that vi+e{0,ipo)) < 0. Consequently the stable 
manifold misses the unstable one for e > small. 

Moreover the stable and unstable manifolds intersect only for a discrete 
set of e, since they vary analytically with e. 

Furthermore equations (fT7|) are reversed by the transformation 

{9,2P)^{-9-n,n-ij). (21) 
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If /3 = 4, the unstable manifold through (— vr, 0) is mapped onto the stable 
manifold through (0,7?) . Hence the stable manifold intersects the line ^ = 
at some point {0,ipo) such that vi+e{0,'ipo)) < and the stable manifold 
misses the unstable one for e > small. 

Moreover the stable and unstable manifolds intersect only for a discrete 
set of e, since they vary analytically with e. 

Similar arguments can be applied to the remaining stable and unstable 
manifolds. This concludes the proof of Proposition 4.1. □ 

Proof of Lemma \0. 1\ Observe that C,^ satisfies the equation 

C^{0,e)= f Fp{rjX^v),e)dri, (22) 

J —n 

where F is given by ( |T71) . For e small we can write 

C^ = Co^(^) + <f(^) + 0(e2). (23) 

We also have 



C^{0) =e + vr. 



(24) 



To compute Ci(^)! we can use the Taylor expansion of ( 1221) with respect to e 
and find 



Standard computations show that 



d^^ ^ ^ /5 cos(r/)sin(r/)cos (^Co(^) 



-F^(r/, Co^(r/), e) = '- ^—^ ^ + O (e) (26) 



and that 

^F,{r^XM,e)=0{e). (27) 
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We can now compute Cf (^) as follows, 



^ rO COS (77) sin (77) cos ( Co (^) 



sin (Co(^) 



ffr,!^ (28) 
df]. 



When /5 = 3, 



3 3 cos (r/) sin (r/) cos (r//2 + 7r/2) 
=2 L sin(,/2 + ^/2) 

^ Jcos(i.)sm(i.)+^. (29) 



1 „\\ A 3 



3 cos - 6* sin 



,2 jj V2 ; 4 

and, in particular, for 6' = 0, we have 



TT 



4 de^ ' ' 2 

When /? = 4, 



Ci'(0) = 7^ = ^C'(0,^). (30) 



M/m 3 r cos r/ sm 7/ cos r/ + 7r 
2 sm [r] + tt) 

= cos (6) sin (6*) + 6* + TT 



and, in particular, for 6 = — 7r/2, we have 



(31) 



This concludes the proof of Lemma 4.1. □ 



V. Escape and Capture Solutions for h = 

We will further study escape (capture) solutions, i.e. the ones for which 
r 00 when t 00 {t —00). From the energy relation we can see 
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that for /i < 0, the radial coordinate r is bounded by the zero velocity curve 
(m = 0,f = 0), so escapes exist only for h > 0. We restrict our analysis to 
the case /i = 0, in which the energy relation ([7]) takes the form 

+ = 2'-"-' + (33) 
With the transformations 

p = r ; V = p 2 v; u = p ^ u, 
the energy relation becomes 

u^ + v^ = 2 + ^p^-\ (34) 

We define the infinity manifold Iq as 

Io = {{p,v,e,u)\ p = 0, Il2+u2 = 2}. (35) 

Since the variable 6* G S*^, the infinity manifold Iq is a torus. 

Remark 0.4. The infinity manifold Iq is independent of the parameter f]. 

After rescaling the time r with the transformation dr = p^ds, the 
equations ([6]) take the form 



dp 
ds 



-pv, 



ds 2 ^ 

- = u, (36) 

1_ 6/3(/i- l)sin2^ 

-;- = H flZ^ 

2 2A^ 



Equations ([31]) and fl5B]) are well defined on the boundary p = 0. Con- 
sequently, the phase space of the coordinates (p, v, 9, u) can be analytically 
extended to contain the manifold Jq. Since dp/ds = for p = 0, this manifold 
is invariant under the flow. 
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Proposition 0.3. All the equilibrium solutions of the flow given by lie 
on the infinity manifold Iq and they form two circles of equilibria given by 

p = o, v = ±V2, ees^ u = 0. 

Proof. It is obvious tliat any point of tlie above circles is an equilibrium 
orbit. If U = in (136|) . by the third equation we have that u = 0, but this 
is a contradiction with the energy relation given by ( |34l) . This proves the 
result. □ 



On the infinity manifold Iq, the equations of motion take the form 

1_2 -, _2 /o 



d9 
ds 

du 1 

— = — uv. 
ds 2 



u, (37) 



We can now prove the following properties. 

Proposition 0.4. The flow on Iq is gradient-like with respect to the v- 
coordinate. 

Proof. From the first equation in (ITTI) . we obtain that > except at 
equilibria, which proves the gradient-like property. □ 

In agreement with the sign of U, and by similarity with the collision 
manifold for /i = 1 studied in Section 4, we also denote the equilibria on Iq 
as C^, respectively. 

Theorem 0.3. On the infinity manifold Iq, the two circles of equilibria 
and C~ are normally hyperbolic. corresponds to a sink, whereas C~ 

corresponds to a source. The escape orbits are the ones having as an 
uj-limit, whereas capture orbits are the ones having C~ as an a-limit. 

Proof. The proof is similar to the one of Theorem 3.1, so we will only sketch 
the main steps. From equation (!34l) . we define 

2b 

G{p, v^e,u)=lf + v^-2- -^P^"' = 0. 
16 



Then I^o is the 3-dimensional manifold given by G~^(0). To study the tangent 
spaces to this manifold at the equilibria, we use as a basis the same vectors 
f fT2|) . Then the linear representation of the vector field fl36|) at any equilibrium 
C+ and C~ is given by the matrix 

















1 








2 



Notice that for C"*", Vq = \f2 and for C~, Vq = —\pl. In the former case 
two eigenvalues are negative and one is zero, whereas in the latter case two 
eigenvalues are positive and one is zero. This completes the proof of the 
normal hyperbolicty and shows the existence of infinitely many escape orbits 
and capture orbits. 

□ 

The next result , which is a direct consequence of Theorem 10. 3^ character- 
izes the flow on the infinity manifold. 

Corollary 0.2. For h = 0, the infinity manifold Iq is foliated by heteroclinic 
orbits between C~ to C"*". 

The flow on Jq given by equations (1371) is easy to draw. Because on Jq, 
u^ + v'^ = 2, we can introduce the angular variable ip with the transformation 

u = V2cosil), V = V^sinip. 

On Jq, the equations of motion take the form 

tp = — 2v^cos^/', 6 = v^cos'0. 

From here we obtain that 

de 

On Iq we can also study the projection of the flow on the v — 6 plane, which 
is given by 

dv _ V2-v^ 
de ~ 2 ' 

whose solution is v{9) = -\/2 sin where is a constant determined by the 
initial condition. 
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VI. Integrability for (3 = 2 

We will further study the problem for (3 = 2 and deal with a Hamiltonian 
([1]) of the form 



With the notation e = /i — 1, the Hamiltonian expressed in polar coordinates 
becomes 

rr ^ /oq^ 

2 2 2r2 r r2(l + e cos2(^)) ' ^ ' 

The corresponding system is integrable since it admits another first inte- 
gral independent of the Hamiltonian, namely 



Indeed, 



and since 



dH.dG dH.dG 



dH2 e sin(2e) dH2 



and 



de r2(l + ecos2(^))2' 
dG esin(2^) (9G 



2 



(42) 



(43) 



de (l + ecos2(^))2' 

the Poisson's bracket is G} = 0, so G and are linearly independent. 

The existence of the integral G is not surprising. Indeed it is well known 
that (see Ref. 19) given an Hamiltonian 

H = ^^ + §r. + UM^ (44) 

the corresponding Hamilton-Jacobi equation 

ds -[ as ds\ ^ , 
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(where S = S{r, 6, t) is the action expressed as function of the coordinates 
and time, dS/dr = pr and dS/dO = pq) can be solved by separation of 
variables if the potential energy is of the form 

b{0) 



U = a{r) + 



(46) 



Since the Hamiltonian is time independent we take S{r, 6, t) = So{r, 9) — Et 
(where ii^ is a constant), and the Hamilton- Jacobi equation for So becomes 



1 /asy 

2 V 9^ J 



a(r) + 



2r2 



(f) -K.) 



E. 



Looking for a solution of the form 



So = Siir) + S2{e), 
we find for and 5'2 the equations 



(47) 



(48) 



1 f dSi 
2 



2h{e) 



2G, 
2G 



(49) 



+ a(r) + ^ = E, 



which define two independent integrals. Solving these equations leads to a 
solution of the Hamilton- Jacobi equations and thus to a general solution of 
the equations of motion. The above technique applies to the Hamiltonian 
H2 and to the additional integral G. 

This approach also shows that the Hamiltonian system given by H2 is 
integrable by quadratures. The same conclusion can be reached by directly 
applying the Liouville-Arnold theorem.^ 



VII. The Collision Manifold for /3 = 2 

For (3 = 2 the equation of motion in McGehee coordinates can be written 

as 



r =uv 

v' =2r'^h + r 

6' =u 

u' =ebsm{2e)A-\ 



(50) 
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where the prime denotes differentiation with respect to r and A = 1 + 
ecos^(^). In McGehee coordinates, the energy relation talces the form 

+ - 2r - 2bA^^ = 2r^h, (51) 

where h is the energy constant. The first integral G can be written as 

g = ^{u^-2bA-'), (52) 

where g is also constant along orbits. 

The vector field fl50|l is analytic on the boundary r = 0, since r no longer 
occurs in the denominators of the vector field. The collision manifold reduces 
to 

C = {{r,e,v,u) -.r = 0, + = 26A"^}. (53) 

This shows that C is homeomorphic to a torus. The restriction of equations 
f l50|l to C yields the system 

v' =0 

e' =u (54) 
u' =ebA~^sm29. 



All nonequilibrium orbits on C are periodic. Comparing the collision 
manifold and vector field above with the corresponding ones in Ref. 6, we 
see that the collision manifold and the flow for (3 = 2 are identical to the 
ones of the anisotropic Manev problem. 



VIII. Heteroclinic orbits for (3 = 2 and h = 

The main goal of this section is to study the infinity manifold for h = and 
the heteroclinic orbits that connect the collision and infinity manifold. First 
notice that for h < the motion is bounded and therefore there is no infinity 
manifold. More precisely, we have the following result. 

Proposition 0.5. If h < the motion is bounded by the zero velocity curve 



ro = ,v = 0,^,M = ) . (55) 
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Proof. Obviously u = and = if and only if + f ^ = 0. Also v = 
implies r' = 0. Using the energy relation we can draw the conclusion that 
+ f ^ = if and only if 2r^/i + 2r + 26A~^ = 0. This quadratic equation 
has the solutions 

-1± Vl-4/i6A-i 

• (56) 



2h 

Since r > and h < the only valid solution is the one with the minus 
sign. This shows that fl55l) is the zero velocity curve. The fact that the 
motion is bounded by this curve follows from the fact that if r > ro then 
u"^ + < 0. □ 

To describe the behavior of the solution at infinity we need to study the 
equations (!36l) with (5 = 2 and = 0, that is 

f) = - pv 
1 „ 

V = W + 1 

e=u ^ ' 

u = — -V u + ehp sin 26'A~^, 

where the dot denotes differentiation with respect to the new time variable 
s. The energy relation is 

u^ + v'^ -2- 2bpA-^ = (58) 

and the other first integral can be written as 

- 2bpA~^ = 2pg. (59) 

The infinity manifold is the set 

Io = {ip,v,e,u)\p = 0, u^ + v^ = 2, u^ = 0}, (60) 

i.e. the points in phase space that satisfy the condition p = 0, the energy re- 
lation and the additional conservation relation, /q is a disconnected manifold 
formed by the union of two disjoint circles of fixed points lying in parallel 
planes. 

The first two equations of the system fl57|) are independent from the 
others, and we would like to determine p and v. If we set v = ±V^, then 
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p = exp(=i=-\/2(s — Sq)) is a solution of the two equations mentioned above. If 
V = ±-\/2 the energy integral ( l58l) gives the condition 



u 



±^/2bpA~\ (61) 



Moreover the previous condition and equation impose g = 0. Differen- 
tiating fl6Tl) with respect to s we obtain 

and using the first equation in f l571) and equation ( 1611) it follows that 

il= -^vu + ebpsm2eA-^. (63) 

This shows that system (157|) admits solutions with U = ±-\/2. 

The other solutions of the first two equations of the system 0571) can be 
found by dividing the second equation by the first. This leads to the equation 

V Iv 1 , 

- = 7; =, 64 

p 2 p pv 

which can be solved by separating the variables. This leads to 



and consequently yields 



^0 



P={J^]iv'-2), (66) 



where po and vq are initial conditions and p > 0, since p < has no physical 
meaning. Now we can prove the following result relative to the heteroclinic 
solutions connecting the infinity and the collision manifolds. 

Theorem 0.4. The solutions whose uj-limit set belongs to the infinity man- 
ifold have the a-limit set contained in the collision manifold. In particular, 
the following properties take place: 
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1. Ifv = \pl (v = —\f2), the above solutions belong to the unstable (stable) 
manifold of one of the periodic orbits on the equator of the collision 
manifold. 



2. IfO< ^/TJk < V2b and ^Jlfk ^ ^26//i with k = po/{vl - 2), the 
above solutions belong to the unstable (stable) manifold of the periodic 
orbits on the collision manifold with v = ^/l/k (v = —^f\Jk). 

3. If y^TJk = a/26/ /i (—^yi/k = ^y2b/^), then the above solutions belong 
to the unstable (stable) manifold of one of the fixed points Aq , (Aq , 



4- If a/I/A; = (-y^l/k = V^), then the above solutions belong to the 
unstable (stable) manifold of one of the fixed points A'^^^^, A^^^ (A~_^i^, 

Proof. To prove (1) observe that v = U/p^/^ = ±-\/2/p^^^- Thus hm^^oo^' = 
hm^^oo±V2/p^^^ = hm^^oo ±V2/p^^^ = hm^^oo^' = 0, since hm^^ooT = 
oo. To prove (2), (3) and (4) we consider the hmit limr^ooV, which with the 
help of equation (!66l) becomes 



V r ^ . r /P + 2fc _i/2 , /I 
hm V = hm —pr = ± hm \ — ; p ' = ±\ -. 

T^OO T— »00 p-*-'^ T—fOO \ k \ k 

Moreover from the energy relation + — 2 — 26p/A = and the fact that 
= (p + 2k) /k it easy to see that 



2b 

k A 



and since p > and > 0, we have 



1 2b 
-< — <2b. 
k - A - 



Consequently we have shown that 



< W ^ < V2b. (67) 



In particular it is clear that when \/TJk = ^J2h[p, {—^Jl/k = a/26/ p) the 
solutions lie on the unstable (stable) manifold of one of the fixed points Aq , 
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{Aq, A~). Similarly when ./TJk = \f% {-\f\Jh = V2b), the solutions 
lie on the unstable (stable) manifold of one of the fixed points ^1^/2' 
(y4~^^2) ^^72) • remaining cases the solutions lie on the unstable (stable) 

manifold of the periodic orbits on the collision manifold with v = \f^^k 

-07fc). □ 



IX. A Perturbative Approach 

In this and the next section we regard the problem in a different way by 
considering the Hamiltonian 

7^, = ip^ - i - A + eh^-^ ^H' + bWiir, 6) + ebW',{r, 6), (68) 

where /3>3/2, 6<^l,e^l. This is the original Hamiltonian H13 (defined 
in ([1])) truncated to the second order. Consider, as in Ref. 5, the parabolic 
solutions of the unperturbed (e = 0) problem (1^ . defined by the Hamilto- 
nian TY" of the classical Kepler problem that are on the zero-energy manifold 
and play the role of homoclinic solutions corresponding to the critical point 
at infinity, i.e., r = 00, f = 0. These solutions satisfy the equations: 

. . k 

r = ± , ^ = (69) 

where 7^ is the (constant) angular momentum and the sign — (resp. +) 
holds for t < (resp. t > 0). From fl69|) we get 

k"^ + r 



We denote by 



±t = -\/2r — fc^ -|- const. 

u = ±2 arctan = h const. 



R = R{t) and 9 = e(t) (71) 



the expressions giving the dependence of r and of 6 on the time t. These 
are obtained by "inverting" the equations flTOj) with the conditions -R(O) = 
"^min = A;^/2 and 6(0) = vr. Let us emphasize that, as in Ref. 5, it is not 
necessary to have the explicit form of these functions. But it is important to 
remark that R{t) is even and G(t) is odd, both in the time variable. 
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The parabolic orbits can also be described in parametric form.^^ If p = 

A;2 ^ 0, we can write 

r = |(l + r/2), t=^r/(l + ^), r/ = tan-. (72) 



We also have 



cos2^ = 2)- (73) 



(l + r/2)2 

We will further use these remarks to apply the Melnikov method. 
X. The Melnikov Method 



Consider the problem defined in ( 1681) . The homoclinic manifold, i.e., the 
set of solution of the unperturbed equation which are doubly asymptotic to 
r = cxD, f = 0, is given for each value k ^ of the angular momentum by the 
2-dimensional manifold described by the family of solutions r = R{t—tQ),'d = 
Q{t — to) + 6*0, where R(t) and 6(t) have been defined in (ITTI) . with arbitrary 
to, ^0- 

It is clear from equation ( l68l) that the first order in b of the perturbation 
(i.e. the term bW^) does not contribute to the Melnikov integrals. Actually 
the perturbed Hamiltonian truncated at the first order is integrable. 

The perturbation resulting from a small anisotropy vanishes as r ^ oo 
since >V^(r, 6') ~ with P > 3/2 satisfies condition (18) in Ref. 5. More- 
over, since the Hamiltonian Tio + is integrable, the terms in 6" for n > 2 
do not contribute either. This allows us to write the first non-vanishing effect 
on the Melnikov integral in the same form as in Ref. 5, with the difference 
that here we can drop the dependence on time, 

Am) ^ r»r^(,) W!)_M±M 

J —oo ^' 



(74) 



dt = 



/•+°° 9wJ(/2(t),e(t) + ^o) , , 

MM = j ^^^^ ^ dt = 0. (75) 

Since the perturbation W| vanishes as t ^ ±oo, the first Melnikov condition 
can be written as 
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Mi(^o) = 




dWl{R{t),Q{t) + eo) 



dt = 0. 



(76) 



The above integral is identically zero because the perturbation is not 
time dependent. This simplifies our discussion since we must only find the 
zeroes of (175|) . 

It is significant to remark, and easy to verify, that these conditions can 
be written also in terms of the first integrals of the unperturbed problem as 



where (...) represents the homoclinic orbit. 

This resembles some properties obtained for the Gylden problem (see 
Refs. 5,16) and is related to the symmetries of the problem. In the Gylden 
problem there is a perturbation that doesn't depend on the angle 6, but 
depends on time. This means that the perturbation destroys the time homo- 
geneity, so the Hamiltonian is not an integral of motion anymore, but doesn't 
destroy the rotational invariance, so the angular momentum is still conserved. 
Therefore only one condition is given by ([71]). However, the anisotropy de- 
stroys the rotational symmetry but not the homogeneity of time, so we are 
left with the condition ( 1751) . 

Here the Melnikov condition for M2 becomes 




(77) 



and 




(78) 




sin[2(e(t) + go)] 
Ritr 



dt = 0. 



(79) 



Using some trigonometry the integral can be written as 



Ii cos 26q + I2 sin 26'o, 



(80) 



where Ii and I2 are defined by 




(81) 
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Recall that R{t) is an even function of time and B(t) is an odd function. 
This implies that the integrand of Ji is an odd function. Therefore Ji = 0, 
and M2 can be rewritten as 



M2(^o) = /2 sin 2^0 



(82) 



Using the parametric form of the parabolic orbits defined in equations 



2]), and observing that dt 



2/3-1^3/2-/3 



2 

+00 



1 + 7]'^)dr], we can write 



00 



2\2 



1 dr]. 



(83) 



Computing the integral, we find that 

2/3-1^3/2-/3^ 



2r(/?-i) 
, r(/5 + i)-r(/3 



r(^-^) 



2(/5-l)/? 



- 1 



{f3 - l)f3 



(84) 



where T{z) is Euler's Gamma function. Thus hiP) is an analytic function 
in /? for (3 > 3/2, since T{z) is analytic for z > 0. Recall that the Gamma 
function can be expressed as 



T{z) = lim 



n\n 



n^oo Z[Z + 1 , 



{z + n) 



(85) 



if 2 7^ 0, —1, —2, —3, . . . Using this form of the Gamma function, and letting 

A = 2/3-2^3/2-/3, we find that 



A lim n 2 x/vr 



2 --—-^"^p_3^^p_i2)...iP+'. + n) 



iP^ - 5/? + 6) 



Av/?r(/5- 



(86) 



(/?-i)(/3-|)(/5-i)r(/3-i) 



which is zero if and only if — 5/3 + 6 = 0, since the Gamma function T{z) 
is always positive for 2 > and therefore the first factor never vanishes. 
Consequently I2 vanishes if and only if — 5/3 + 6 = 0, namely for /3 = 2 or 
P = 3, see Figure [3l 

We have thus proved the following result. 
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Figure 3: I2/A as a function of p. 



Theorem 0.5. For every (3 > 3/2, 7^ 2,3, the system given by the Hamil- 
tonian ( fg^) exhibits chaotic dynamics on the zero energy manifold. 

This type of chaotic behavior is induced by a chain of infinitely many 
intersections of the positively and negatively asymptotic sets to the critical 
point at infinity. The Smale-Birkhoff theorem does not directly apply to 
this situation, which is degenerate. But it is well known that the existence 
of Smale horseshoes and positive topological entropy can arise in the case 
of nonhyperbolic equilibria (see for example Ref. 4, where the problem is 
studied in the case of area-preserving diffeomorphisms) . 
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